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The Chapman-Kolmogorov equation with fractional integrals is derived. An integral of frac- 
tional order is considered as an approximation of the integral on fractal. Fractional integrals can be 
used to describe the fractal media. Using fractional integrals, the fractional generalization of the 
Chapman-Kolmogorov equation is obtained. From the fractional Chapman-Kolmogorov equation, 
the Fokker-Planck equation is derived. 

1 Introduction 

Integrals and derivatives of the fractional order goes back to Leibniz, Liouville, Riemann, 
Grunwald, and Letnikov [T|. Fractional analysis has found many applications in recent 
studies in mechanics and physics. The interest in fractional integrals and derivatives has 
been growing continually during the last few years because of numerous applications. In a 
fairly short period of time, the list of such applications has become long. It includes chaotic 
dynamics mechanics of fractal and complex media [H El E], physical kinetics [21 H E], 
plasma physics [HI HOI IH] , astrophysics [12] , long-range dissipation [131 [E] , non-Hamiltonian 
mechanics [151 IIH] > cind long-range interaction [T71 [HI [19] . 

The natural questions arises: What could be the physical meaning of the fractional inte- 
gration? This physical meaning can be following: the fractional integration can be considered 
as an integration in some noninteger-dimensional space. If we use the well-known formulas 
for dimensional regularizations [20], then we get that the fractional integration can be con- 
sidered as an integration in the fractional dimension space [TB] up to the numerical factor 
r(a/2)/[2W2r(a)]. This interpretation was suggested in Ref. [15]. Fractional integrals can 
be considered as approximations of integrals on fractals [211 [22]- In Ref- [22], authors proved 
that integrals on a net of fractals can be approximated by fractional integrals. Using frac- 
tional integrals, we derive the fractional generalization of Chapman-Kolmogorov equation 
[231 [21] • In this paper, the generalization of the Fokker-Planck equation for fractal media is 
derived from the fractional Chapman-Kolmogorov equation. 

In Sec. 2, a brief review of the Hausdorff measure, the Hausdorff dimension and integra- 
tion on fractals is carried out to fix notation and provide a convenient reference. The connec- 
tion of integration on fractals and fractional integration is discussed. We derive the fractional 
generalization of the average values equation. In Sec. 3, the fractional Chapman-Kolmogorov 
equation is derived by using fractional integration. In Sec. 4, the fractional Fokker-Planck 
equation for the fractal media is derived from the suggested fractional Chapman-Kolmogorov 
equation. The stationary solutions of the Fokker-Planck equation for fractal media are de- 
rived. 
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2 Integration on Fractal and Fractional Integration 



Fractals are measurable metric sets with a non-integer Hausdorff dimension. Let us consider 
a brief review of the Hausdorff measure and the Hausdorff dimension in order to fix notation 
and provide a convenient reference. 



2.1 Hausdorff measure and Hausdorff dimension 

Consider a measurable metric set {W, fin)- The elements of W are denoted hj x,y, z, and 
represented by n-tuples of real numbers x = {xi, X2, Xn) such that W is embedded in R^. 
The set W is restricted by the conditions: (1) W is closed; (2) W is unbounded; (3) W is 
regular (homogeneous, uniform) with its points randomly distributed. 

The metric d{x, y) as a function of two points x and y (zW can be defined by 



Xi 



i=l 

The diameter of a subset E (ZW <Z is 

d{E) = diam{E) = sup{d{x,y) : x,y E E}. 

Let us consider a set {Ei} of non-empty subsets Ei such that dim{Ei) < e, Wi, and 
W C U^=i Ei. Then, we define 

aEi,D) = uj{D)[dtam{E,)f = uj{D)[d{Ei)f . (2) 

The factor uj{D) depends on the geometry of Ei, used for covering W. If {Ei} is the set of 
all (closed or open) balls in W, then 

The Hausdorff dimension D of a subset E cW is defined [25], [26] by 

D = dimniE) = snp{d G R : fJ^niE, d) = oo} = in{{d G R : fJ^niE, d) = 0}. (4) 

From (jl]), we obtain fiH^E, d) = for d > D; and finiE, d) = oo for d < D. 
The Hausdorff measure of a subset E (ZW is [251 126] : 



lJiH{E,D) = \im mi{Y,i{Ei,D)-. E (z[jE,, d{E,) < e Vz}, (5) 
^^^i^-i i=i i 

or 

oo 

fXHiE,D) =u;iD) Jim iniJ2[d{Er)]''. (6) 
If E C and A > 0, then /i/f(AE, £)) = X^finiE, D), where = {Xx, x G E}. 
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2.2 Function and integrals on fractal 

Let us consider the functions 

oo 

f{x) = J2(3aEA^), (7) 

i=l 

where xe is the characteristic function of E: Xe{x) = 1 if x E E, and Xe{x) = if x ^ E. 
The Lebesgue-Stieltjes integral for ([7]) is defined by 



„ oo 

Jw ^ 



Therefore 



/ f{x)dMx) = y.fi^^)m.D) = uj{D) hm J2fi^^ME^)]''■ (9) 

JW cx(_fc/i)^U g a(_fc/ij— >U g 

It is always possible to divide into parallelepipeds: 

Ei^...i„ = {{xi, ...,Xn) eW : Xj = {ij - l)Axj + aj, < aj < Axj, j = 1, ...,n}. (10) 
Then 



n 



dfiH{x)= \im e(i?,,..,„,D)= hm mAx,)''/- = l[d^/-x,. (11) 

The range of integration W may also be parametrized by polar coordinates with r = d{x, 0) 
and angle Q. Then E^^^ can be thought of as spherically symmetric covering around a center 
at the origin. In the limit, the function ^{Er^n,D) gives 

diiH{r,n)= lim i{Er^,D) = dVt^'h^^^dr. (12) 

Let us consider /(x) that is symmetric with respect to some point xq G W , i.e. /(x) = 
const for all x such that d{x,Xo) = r for arbitrary values of r. Then the transformation 

W : x^x' = x-Xo (13) 

can be performed to shift the center of symmetry. Since W is not a linear space, the 
transformation (fT3l) need not be a map of W onto itself, and (fT3l) is measure preserving. 
Then the integral over a D-dimensional metric space is 

/ fdfiH = KD) / /(r)r^-Mr, (14) 

where X{D) = /V{D /2). This integral is known in the theory of the fractional calculus 
[1]. The right Riemann-Liouville fractional integral is 



oo 



(/^'/)W = ^/ {x-zr-'f{x)dx. (15) 



Then Eq. (fT^ is reproduced by 



/ 



„/*« = ^(^(/-"/)(0). (16) 

Equation f[T^ connects the integral on fractal with integral of fractional order. This result 
permits to apply different tools of the fractional calculus [1] for the fractal medium. As a 
result, the fractional integral can be considered as an integral on fractal up to the numerical 
factor T{D/2)/[2tt''/^T{D)]. 

Note that the interpretation of fractional integration is connected with fractional di- 
mension [12]. This interpretation follows from the well-known formulas for dimensional 
regularization [20] . The fractional integral can be considered as an integral in the fractional 
dimension space up to the numerical factor T{D/2)/[27i^^^T{D)]. In Ref. [21], it was proved 
that the fractal space-time approach is technically identical to the dimensional regularization. 



2.3 Properties of integrals 

The integral defined in Eq. (Q satisfies the properties: 

(1) Linearity: 

/ (a/i + bf2) dfiH = a fi dfin + b /2 dfin, (17) 
Jw JW JW 

where /i and /2 are arbitrary functions; a and b are arbitrary constants. 

(2) Translational invariance: 

/ f{x + xo)diJ,H{x)= f{x)d^H{x) (18) 

Jw JW 

since dfinix — xq) = dfinix) as a consequence of homogeneity (uniformity). 

(3) Scaling property: 

/ /(Ax)rf/iH(x) = A-^/ fix)dfiHix) (19) 

JW JW 

since d^H^x/X) = X^^dfinix). 

It is well-known [271 120] that conditions f|T7j) - (fT9l) define the integral of the function 
f{x) = exp(— ax^ -|- bx) up to normalization: 

/ exp(-aa;2 + bx) c//i//(x) = vr^/^^"^/^ exp(6V4a). (20) 

JW 

For 6 = 0, Eq. fl20|) is identical to the result that can be derived from f|T6|) and is obtained 
directly without conditions f|T7j) - f|T9l) . 
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2.4 Fractional average values 

The usual average value 

/ + 00 
A{x)p{x)dx (21) 
-oo 

can be written as , _ 

A{x)p{x)dx + / A{x)p{x)dx. (22) 

Using 

the average value fl22]) can be present by 

<A>i=(J^Ap)(y) + (jMp)(y). (25) 
The fractional generalization of Eq. fl25p is 

< A (y) = {II Ap) (y) + (rAp) (y) . (26) 
Equation (126|) can be rewritten as 

<A>^{y)= [{Ap){y-x) + {Ap){y + x)]dp^{x), (27) 



where 

c^/^a = p. ^ = -pT^- 28 
1 (a) ai [a) 

Here, we use 

= /3(x)(x)" = sc/n(x)|x|", (29) 

where j3{x) = {sgn{x))'^~^. The function sgn{x) is equal to +1 for x > 0, and —1 for x < 0. 
To have the symmetric limits of the integral, we consider (127|) in the form 

X r+oo 

<A>^{y) = - [{Ap){y - x) + {Ap){y + x)]dp,(x). (30) 

Z J— oo 

If a = 1, then we have the usual equation for the average value. 

Let us introduce some notations to simplify Eq. ( l30l) . We define the integral operators 

ilfix) = I r^lfix) + f{-x)]dp^ix). (31) 

/ J —oo 

Then fl5Ul) has the form 

< A >a= I^A{x)p{x). (32) 

We will use the initial points that are set to zero {y = 0). Note that the fractional normal- 
ization condition is a special case of this definition of average values: < 1 >a= 1- 
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3 Fractional Chap man- Kolmogorov (FCK) Equation 

The Chapman-Kolmogorov equation [23l [23] may be interpreted as the condition of consis- 
tency of distribution functions of different orders. Kolmogorov [231 [2l] derived a kinetic equa- 
tion using a special scheme and conditions that are important for kinetics. Let W{x, t; xq, to) 
be a probability density of having a particle at the position x at time t if the particle was at 
the position xq at time to < t. 

Denote by p{x, t) the distribution functions for the given time t. Let us consider two 
well-known identities 

/ + 00 /■ + 00 

dx' W{x,t\x',t')p{x',t'), / p{x,t) = l. (33) 
-oo J—oo 

Using the notation fl3Tl) . we can rewrite fl33l) in the form 

p{x, t) = II, W{x, t\x\ t')p{x', t'), II p{x, t) = 1. 
Then the fractional generalization of fl33|) is 

p{x, t) = I^, W{x, t\x', t')p{x', f). (34) 

This equation is the definition of conditional distribution function W{x,t\x\t') referring to 
different time instants. The normalization conditions for the functions W{x^t\x' ,t') and 
p(x, t) are 

1° W{x, t\x\ t') = 1, 4" p(x, t) = 1. (35) 

Substituting into the right-hand side of Eq. (IMll the value of p{x',t') expressed via the 
distribution p(xo,to) at an earlier time, 

p{x',t') = W{x',t'\xo,to)p{xo,to), (36) 
we obtain the integral relation which includes the intermediate point x', 

p{x,t) = i:, W{x,t\x',t')W{x',t'\xo,to)p{xo,to). (37) 
Using Eq. fl371) . and Eq. (134]) in the form 

p{x, t) = W{x, t|xo, to)p(a;o, to), (38) 
we derive a closed equation for transition probabilities 

/° W{x,t\xo,to)p{xo,to) = i:. 1° W{x,t\x\t')W{x' ,t'\xo,to)p{xo,to)■ 
'&m.ce the equation holds for arbitrary p(a;o,to)? we may equate the integrand. As the result, 
we obtain the fractional Chapman-Kolmogorov (FCK) equation 

W{x, t\xo, to) = 1°, W{x, t\x', t')W{x', t'\xo, to). (39) 

This equation can be used to describes the Markov-type process in the fractal medium that 
is described by the continuous medium model [5]. 
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4 Fokker-Planck Equation from FCK Equation 

4.1 Derivations of Fokker-Planck equation 

Let us consider the fractional average value fl32|) . Using Eq. flM|) in the form 

t) = W{x, t\xo, to)p{xo, to), (40) 

we get 

< A >,= 4" A{x) 4"^ iy(x,t|xo,to)p(a;o,to). (41) 
We can rewrite Eq. (HT]) as 

< A >,= 4"^ p(xo, to) 4" ^(a;) W^(a;, t|xo, to)- (42) 
We assume that A = yl(x"), and use the Taylor expansion 

Aix-) = A(x^ + Ax") = Aix-) + (^^)^ Ax" + I (^^)^ (Aa:")^ + (43) 

where x" = S(?n(x) |x|" is defined by Eq. ( !29l) . Ax" = x" — Xg, and 

9 |xP~" d 



dx"' a dx 



(44) 



If we use the usual Taylor expansion, then the integration by parts in Eq. (1421) is more 
complicated. For the expansion (H3l) . the integration by parts in (l42l) can be realized in the 
simple form, 



dx"" J-oo aT(a) dx 



a 



-oo 



{B{x)A{x)tZ- I ^A(x"^^^^^) 



oo 



aT(a) dx" 



Substituting Eq. (HSj) in Eq. (Hg), we get 

< A >„= 4"^ A(x^)p(xo,to)4" W{x,t\xo,to) + 

+Co f^4^) p(xo,to)4" Ax"iy(x,t|xo,to) + 



dx'^ 



/ Xo 



+ ^4°o (^^) p(xo,to)4"(Aa;")2W^(x,t|xo,to) + ... (45) 
ions: 

P„(xo,t,to) = i: (Ax")'^W^(x,t|xo,to). (46) 



/ XQ 

Let us introduce the functions: 
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Using (H6il and (]35ll . equation (H5l) gives 

<^>«=4"o ^K)pK,to) + 4"o (^4^1 p(a:o,to)A(xo,t,to) + 



+ 2^=^n"(9^J P(^o,to)P2(xo,t,to) + ... . (47) 
Substitution of (!32|) in the form 

< A>„=/° A(x^)p(xo,t), 

into Eq. (H71) gives 

4" ^(2^0) (p(a;o,t) - p(a;o,to)) = 

= C P(^o,to)Pi(^o,t,to) + ^4" (1^^) p(a;o,to)P2(xo,t,to) + ... . (48) 

Then we use so-called Kolmogorov condition [23l [2l], and assume that the following finite 
limits exist: 

Pl{x,t,to) / , \ r P2{x,t,to) , / , X Pn{x,t,to) 

hm = aixAo)., lim = MxAn), hm = U, 

where n = 3, 4, and At = t — to- It is due to the Kolmogorov conditions that irreversibility 
appears at the final equation. Multiplying both sides of Eq. (HHl) by 1/At and consider the 
limit At — i> 0, we obtain 

/ dp{xo,t)\ 



4;o ^(^0/ 



\ dt 



to 



4" (^^j p(xo,to)a(a:o,to) + -4°, [-Jd^J 



Integrating by parts, we obtain 



{dx 



p{x,t)b{x,t) = 4" A{x' 



d'^{p{x,t)b{x,t)) 



{dx 



a\2 



(50) 



Here, we use 



Then 



lim p(x, t) = 0. 

a;— > it 00 



Ja A( a. ( dpix,t) d{p{x,t)a{x)) 1 d^pjx, t)b{x)) \ ^ 

" ^ '[ dt ^ dx- 2 {dx-y i 
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Since the function A = 74(x") is an arbitrary function, then we have 

dp{x,t) d{p{x,t)a{x,t)) l d'^{p{x,t)b{x,t)) _ 

dt dx^ 2 {dx^f ' ^ ^ 

that is the Fokker-Planck equation that corresponds to the FCK equation. This equation is 
derived from the fractional generalizations of the average value and fractional normalization 
condition, which use the fractional integrals. 

4.2 Stationary solutions 

For the stationary case, the Fokker-Planck equation flSip is 

d{p{x, t)a{x, t)) 1 d^{p{x, t)b{x, t)) 



2 ((9x")2 

This equation can be rewritten as 



0. (52) 



d f , , , , 1 d(p(xA)b(x,t))\ , , 

' p(x, t)a(x, t) - - \[} " = 0. (53) 



ax° V ' ' ' ' 2 ax° 

Then 



p(..i)a(..t)-i ^<''™-") ^c.,,.t. (54) 



Supposing that the constant is equal to zero, we get 

d{p{x,t)b{x,t)) 2a(x, t) 



(p(x,t)6(x,t)), (55) 



ax" 6(x, t) 

The solution of (l55l) is 

ln(p(x, t)b(x, t)) = f ^^^^dx" + const. (56) 
J b{x,t) 

As the result, we obtain 



N f a(x, t) , „ , , 



b{x,t) J b{x,t) 

where the coefficient N is defined by the normalization condition. 

Let us consider the special cases of the solution (IFH) . 
(1) If a(x) = k and b{x) = —D, then the Fokker-Planck equation fl5T|) has the form 

ap(x,t) ap(x,t) D d^p{x,t) _ 
dt ^ dx- ^2 (ax")2 ^""^ 
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and the stationary solution is 



p{x,t) = Niexp ( -^^1^ ) • (59) 



D 



(2) If a{x) = k\xf and b{x) = -D, then 



dpix, t) a^lVM) , D d^p{x,t) _ 



at 2 (9x" 

The stationary solution is 

2ak\x 



If a + /3 = 2, we have 



(3) If 



and b = —D, then 



pi,,t)^N,_e.p{-=^^j. (61) 



p(x, t) = N^exp -^x^ 1 • (62) 



, , dUix) \x\^-"dU{x) 

aix) = — = — , 

ox" a ox 



p{x,t) = N^exp ( 



Let us consider Eq. flHT]) with a{x) = k\x\°' and b = —D. The general solution can be 
presented as 



+ 00 



where 



The stationary solution is 

1/2 
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5 Conclusion 



The concept of fractional integration provides an approach to describe the fractal media. The 
fractional integrals can be used in order to formulate the dynamical equations in the fractal 
media. The fractional integration approach is potentially more useful for the physics of fractal 
media than traditional methods that use the integer integration. Using fractional integrals, 
we derive the fractional generalization of the Chapman-Kolmogorov equation. This equation 
can be used to describes the Markov-type process in the fractal medium that is described by 
the continuous medium model [5]. The fractional Chapman-Kolmogorov equation can have 
a wide application since it uses a relatively small number of parameters that can define a 
fractal medium of great complexity and rich structure. In this paper, we derive the Fokker- 
Planck equations for the fractal media from the suggested fractional Chapman-Kolmogorov 
equation. 
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